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ABSTRACT 


Since  the  1960s,  comparing  a  Rayleigh-wave  magnitude,  M,  to  the  body-wave  magnitude,  mb  has  been  a  robust  tool 
for  the  discrimination  of  earthquakes  and  explosions  (e.g.,  Ms:mb).  In  this  article,  we  apply  a  Rayleigh-wave  formula 
as  is  to  Love  waves  and  examine  the  possibilities  for  discrimination  using  only  surface  wave  magnitudes  (e.g., 
Mj(Rayleigh):Ms(Love)).  To  calculate  the  magnitudes  we  apply  the  time-domain  magnitude  technique  called 
MS(VMAX)  developed  by  Russell  (2006)  to  Rayleigh  and  Love  waves  from  explosions  and  earthquakes.  Our  results 
indicate  that  for  the  majority  of  the  earthquakes  studied  (>75%),  the  M,(VMAX)  obtained  from  Love  waves  is 
greater  than  the  estimate  from  Rayleigh  waves.  Conversely,  79  of  82  nuclear  explosions  analyzed  (96%)  have 
network-averaged  MS(VMAX)-Rayleigh  equal  to  or  greater  than  the  M,(VMAX)-Love. 

We  use  logistic  regression  to  develop  an  Mj(Rayleigh):Ms(Love)  discriminant.  Cross-validation  analysis  of  the  new 
discriminant  correctly  identifies  57  of  82  explosions  and  246  of  264  earthquakes  while  misidentifying  22  explosions 
as  earthquakes  and  1 1  earthquakes  as  explosions.  The  majority  of  these  misidentified  earthquakes  are  either  deep 
(sub-crustal)  earthquakes  or  the  events  with  normal  or  thrust  focal  mechanisms.  Further  comparative  research  is 
planned  for  M,(Rayleigh):M,(Love)  versus  Ms:mb  using  common  data. 

We  expect  that  M v ( R a y  I e i g h ) : A7 v ( L o v e )  will  contribute  significantly  to  multivariate  event  identification;  however,  it 
does  not  have  the  same  population  separation  that  has  been  historically  observed  for  the  Ms:mb  discriminant.  Results 
also  suggest  that  incorporation  of  Love  waves  into  the  analysis  requires  a  re-examination  of  the  period  limits 
currently  used  for  the  Af/VMAX)  technique.  While  the  A7S(V MAX)- Rayleigh  method  is  currently  operational  at 
different  data  centers  using  periods  between  8  and  25  seconds,  we  believe  that  future  processing  should  be  extended 
to  40  seconds,  especially  in  regions  with  deep  earthquakes  and  complex  paths. 
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OBJECTIVES 


The  objective  of  this  research  is  to  apply  an  Ms  formula,  originally  developed  and  applied  to  Rayleigh  waves,  to  both 
Love  and  Rayleigh  waves.  During  the  past  year,  we  have  applied  the  method  in  the  same  manner  to  both  phases  for 
three  different  earthquake  datasets  as  well  as  a  global  dataset  of  nuclear  explosions.  We  examine  whether  improved 
discrimination  is  possible  by  combining  the  Love  and  Rayleigh  wave  magnitudes.  Finally,  we  discuss  possible 
methods  for  improving  the  analysis  (e.g.,  using  longer  surface  wave  periods,  different  attenuation  corrections,  etc.) 
based  on  results  of  this  study. 

RESEARCH  ACCOMPLISHED 


Love  Wave  Magnitude  Estimation 

We  evaluate  applicability  ofMJVMAX)  (Variable-period,  MAXimum  amplitude  surface  wave  magnitude 
estimation)  to  Love  waves.  The  formula  for  MS(V MAX)  was  developed  by  Russell  (2006)  while  the  measurement 
technique,  which  is  currently  in  use  at  the  United  States  Geological  Survey  as  Ms_VX,  was  developed  by  Bonner  et 
al.  (2006).  MJVMAX)  was  developed  for  Rayleigh  waves  measured  at  variable  periods  between  8  and  25  seconds. 
It  is  defined  as  follows: 

MS(VMAX )  =  log(A)  +  ilog {sin  A)  +  0.0031  A  -  log (/c)  -  0.43  -  0.66  log  (^),  (1) 


where  T0  =  20  sec  is  the  reference  period,  A  is  the  great  circle  distance  in  degrees,  fc  is  the  comer  frequency  of  the 
filter,  the  constant  (0.43)  was  obtained  for  zero-phase,  third  order  Butterworth  filter.  The  second  term  of  Equation  1, 

-log(sm  A)  ,is  a  correction  for  the  geometrical  spreading,  the  third  term,  0.0031  A,  is  a  period-dependent 

attenuation  correction  and  the  fifth  term,  —0.66  log  is  a  period-dependent  excitation  correction. 

For  this  paper,  we  apply  Equation  1  as  is  to  Love  waves.  We  use  the  same  processing  for  the  Love  waves  as 
Bonner  et  al.  (2006)  designed  for  Rayleigh  except  that  for  the  Love-wave  magnitude  estimates  we  filter  the 
transverse  components.  Examples  of  Butterworth  filtering  for  Rayleigh  and  Love  waves  are  shown  in  Figure  la,  b. 
The  data  are  filtered  at  center  periods  of  8,  9,  10,... 25  seconds,  the  maximum  amplitude  at  each  period  is  measured, 
and  Equation  1  is  used  to  form  18  different  magnitude  estimates  for  each  station  (Figure  lb).  The  magnitude  at  the 
period  of  the  maximum  amplitude  is  assigned  as  the  MS(V MAX)  for  a  particular  station,  and  combined  with  other 
stations  to  form  a  network  average  for  an  event.  For  this  study,  analysts  (the  first  two  co-authors)  identified  all 
Rayleigh-  and  Love-wave  phases;  however,  we  are  currently  working  on  automated  methods  to  identify  the  phases 
and  measure  the  amplitudes. 

Earthquakes.  Equation  1  was  applied  to  estimate  M,(V MAX)  for  both  Rayleigh  and  Love  waves  for  three  separate 
earthquake  datasets  (Figure  2).  The  first  dataset  (See  Figure  2a)  included  109  events  located  in  the  Middle  East  with 
the  body  wave  magnitudes  ranging  between  3.8  and  6.1.  The  database  samples  a  variety  of  different  focal 
mechanisms.  The  stations  used  to  estimate  surface  wave  magnitudes  are  distributed  throughout  Eurasia  with 
distances  ranging  from  approximately  83  to  over  10000  km.  The  data  for  these  stations  were  obtained  from  the 
Incorporated  Research  Institutions  (IRIS)  in  Seismology  Data  Management  Center  (DMC),  corrected  for  the 
instrument  response  to  displacement  in  nanometers,  and  rotated  to  transverse,  radial,  and  vertical  components.  The 
Love  wave  magnitudes  were  estimated  from  the  transverse  data,  while  the  Rayleigh  wave  estimates  were  obtained 
using  the  vertical  data.  The  results  are  plotted  in  Figure  3a  and  show  that  M,(VMAX)-Love  exceeds  or  is  equal  to 
MS(VMAX)-Rayleigh  for  82  out  of  109  events  (75%).  The  dominant  periods  of  the  measurements  are  approximately 
21-22  sec  for  Rayleigh  waves  and  25  sec  for  Love  waves.  The  interstation  standard  deviation  averaged  0.22 
magnitude  units  (m.u.)  for  both  Rayleigh  and  Love  waves. 

The  second  dataset  included  31  earthquakes  occurring  in  the  Korean  Peninsula  and  surrounding  regions  (Figure  2c). 
These  events  ranged  in  size  between  3.2  <  Mw  <5.1  with  the  focal  mechanisms  being  predominantly  strike-slip 
(http://www.eas.slu.edu/eqc/eqc_mt/MECH.KR,  last  accessed,  April  2011).  The  distances  to  the  three-component 
stations  recording  these  events,  mainly  Korean  Meteorological  Administration  (KMA)  and  some  Global 


437 

Approved  for  public  release;  distribution  is  unlimited. 


2011  Monitoring  Research  Review:  Ground-Based  Nuclear  Explosion  Monitoring  Technologies 


Seismographic  Network  (GSN)  stations,  ranged  from  55  km  to  1900  km.  Similar  to  the  Middle  Eastern  events,  the 
majority  of  these  events  (25  out  of  3 1,  or  80%)  had  Ms(VMAX)-Love  exceeding  or  equal  to  the  M,(VMAX)- 
Rayleigh,  which  is  expected  for  a  strike-slip  mechanism.  The  dominant  periods  of  the  measurements  for  Rayleigh 
waves  was  less  than  13  seconds;  however,  the  Love  wave  magnitudes  were  uniformly  sampled  between  periods  8 
and  20  seconds.  The  interstation  magnitude  standard  deviation  for  the  Rayleigh  and  Love  waves  averaged  0.1 1  and 
0.22  m.u.,  respectively. 


a)  Rayleigh  Waves  b)  Love  Waves 
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c)  Magnitude  Estimation 


Figure  1.  Examples  of  the  Ms  (YMAX)  technique  applied  to  a)  Rayleigh  and  b)  Love  waves  from  a  Nevada 
Test  Site  explosion  (Hoya).  The  Butterworth  filters  are  computed  at  center  periods  between  8  and 
25  seconds  (not  all  filter  panels  are  shown  in  this  figure).  The  maximum  amplitude  in  each  filter 
band  in  a  Rayleigh  and  Love  wave  group  velocity  window  (small  vertical  lines)  is  input  into 
Equation  1  and  18  different  magnitudes  c)  are  estimated.  The  magnitude  at  the  period  of  maximum 
amplitude  (shown  as  a  star)  is  used  as  the  final  Ms  (VMAX)  for  a  station  and  combined  with  others 
for  a  network  average. 
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a)  Middle  Eastern  Earthquakes 


b)  Korean  Peninsula  Earthquakes 


c)  Italian  Earthquakes 


Figure  2.  Map  of  the  seismic  events  for  which  Ms 
(VMAX)-Love  and  Rayleigh  were 
estimated  in  the  a)  Middle  East,  b) 
Korean  Peninsula  region,  and  in  c) 
central  Italy. 


b)  Korean  Peninsula  Earthquakes 
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c)  Italy  Earthquakes 


Figure  3.  Ms(VMAX)-Love  versus  Ms  (VMAX)- 
Rayleigh  for  earthquakes  in  the  a) 
Middle  East,  b)  Korean  Peninsula 
region,  and  in  c)  central  Italy. 


The  third  dataset  focused  on  the  damaging  L'Aquila  earthquake  (6  April  2009  Mw=6.1)  and  its  aftershocks 
(Figure  2c).  We  have  estimated  M,(VMAX)  for  125  Italian  earthquakes  with  2.8  <  Mu,<  6.1  using  Istituto 
Nazionale  Geofisca  e  Vulcanologia  (INGV)  stations  at  distances  ranging  from  50  to  414  km.  We  include  these  data 
in  our  study  because  the  dominant  focal  mechanism  suggests  NW/SE  trending  normal  faults 
(http://www.eas.slu.edu/eqc/eqc_mt/MECH.IT/laquila.png,  last  accessed  April  2011;  Herrmann  and  Malagnini, 
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2009).  Figure  3c  shows  that  for  the  majority  of  these  events  (100  or  80%)  M,(VMAX)-Love  was  on  average  0.2 
m.u.  larger  than  the  M,(VMAX)-Rayleigh,  which  is  unexpected  for  the  dip-slip  focal  mechanisms.  The  dominant 
period  of  the  measurements  was  8  seconds  for  Rayleigh  waves  and  between  8  and  12  seconds  for  the  Love  waves. 
The  interstation  standard  deviation  for  the  Rayleigh  waves  averaged  0.17  m.u.,  which  was  slightly  lower  than  for  the 
Love  waves  (0.20  m.u.). 

The  percentage  of  the  events  with  higher  M,(VMAX)-Love  is  slightly  lower  for  the  Middle  East  dataset  than  for 
Korea  and  Italy  (74%  vs.  80%).  Possible  explanations  include  deeper  events  as  well  as  more  variety  in  the  focal 
mechanisms  for  the  Middle  Eastern  dataset.  The  interstation  standard  deviation  is  slightly  higher  for  the  Middle  East 
data,  which  most  likely  results  from  more  laterally  heterogeneous  structure.  Another  peculiarity  of  the  Middle  East 
dataset  is  the  longer  dominant  periods  at  which  MS(V MAX)  is  calculated  for  both  Rayleigh  and  Love  waves.  Figure 
4  shows  the  histograms  of  the  dominant  periods  for  the  Middle  East,  Korea,  and  Italy.  Excitation  due  to  depth  alone 
cannot  explain  this  period  increase,  because  a  similar  feature  is  observed  for  a  nuclear  explosion  detonated  in  this 
study  region  (EVID  19980528101600).  The  large  number  of  Love  and  Rayleigh-wave  observations  at  25  seconds 
represents  an  edge  effect  associated  with  the  long  period  limit  in  the  current  processing.  Increasing  this  limit  to 
40  seconds  will  be  discussed  later  in  the  paper. 
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Figure  4.  Histograms  of  the  periods  of  maximum  amplitudes  for  Rayleigh  and  Love  waves  in  the  a)  Middle 
East,  b)  Korean  Peninsula  region,  and  in  c)  central  Italy.  For  the  Middle  East  dataset,  varied  focal 
mechanisms,  depths,  and  complex  regional-to-teleseismic  propagation  paths  lead  to  longer  period 
magnitude  estimates.  For  the  Korean  and  Italian  datasets,  the  events  are  shallow  and  have  shorter, 
less  complex  propagation  paths  leading  to  more  short-period  magnitude  estimates. 
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Explosions.  We  have  also  estimated  the  MS(VMAX)  for  Rayleigh  and  Love  waves  from  82  nuclear  explosions  at 
many  different  test  sites  (Figure  5).  Our  working  hypothesis  was  that  the  Love  wave  magnitudes  should  be  smaller 
than  the  Rayleigh  wave  estimates  for  explosions.  This  was  certainly  the  case  for  all  analyzed  events  at  the  Nevada 
National  Security  Site  (NNSS) — formerly  the  Nevada  Test  Site — where  M,(VMAX)-Rayleigh  averaged  0.4  and 
m.u.  larger  than  Ms(VMAX)-Love.  There  were  some  events  with  large  Love  wave  magnitudes  from  the  Shagan  Test 
Site;  however,  the  Rayleigh  wave  magnitudes  were  on  average  0.21  m.u.  larger  than  the  Love  wave  estimates.  For 
7  Lop  Nor  explosions  the  MS(VMAX)-Rayleigh  are  slightly  larger  than  Mf(VMAX)-Love,  except  for  one  anomalous 
event  (EVID  19920521045947  in  ES  Table  3)  which  had  a  Love  wave  magnitude  0.27  m.u.  larger  than  the  Rayleigh 
magnitude.  Our  dataset  also  included  the  1998  Pakistan  nuclear  test,  which  had  a  Love  wave  magnitude  slightly 
larger  than  the  Rayleigh  magnitude,  and  the  1 998  Indian  nuclear  explosion,  which  had  a  larger  Rayleigh  magnitude 
(by  ~0.2  m.u.). 


We  were  unable  to  measure  Love  waves  using  openly  available  data  for  the  2006  North  Korean  nuclear  explosion 
(MS(VMAX)-Rayleigh=2.9).  Based  on  background  noise  levels,  we  conclude  the  Mf(VMAX)-Love  must  have  been 
less  than  2.5  (similar  conclusion  reached  by  Kohl  et  ah,  2011).  For  the  2009  event,  the  Rayleigh  MS(VMAX)=3.7 
exceeded  the  Love  MS(V MAX)  by  0.5  m.u.  As  mentioned  previously,  the  Korean  events  had  large  Rayleigh  Ms 
estimates  compared  to  mb,  however,  the  Love  waves  magnitudes  are  much  smaller  and  provide  added  discrimination 
information. 


Explosions 


Figure  5.  Ms(VMAX)-Love  versus  MS(VMAX)-Rayleigh  for  nuclear  explosions. 

Event  Identification  with  Logistic  Regression 

The  observed  differences  in  the  Love  and  Rayleigh  wave  magnitudes  between  earthquakes  and  explosions  led  us  to 
an  idea  that  a  surface  wave  discriminant  could  be  developed  without  incorporation  of  an  mb.  For  regional  events, 
mb(Pn)  is  often  difficult  to  determine,  and  there  may  be  geophysical  structural,  data  center  measurement 
(e.g.,  Murphy  et  ah,  1997),  and  data  censoring  biases  that  complicate  the  Ms\mb  interpretation.  We  decided  to  test  for 
a  possible  MS\MS  discriminant  using  logistic  regression  (Press  and  Wilson,  1978). 

Logistic  regression  models  the  conditional  probability  that  an  event  is  an  explosion  given  a  regression  function  of 
event  magnitudes  X .  The  calibrated  model  gives  the  best  linear  combination  (regression  model)  of  magnitudes 
(the  discriminant)  that  best  agrees  with  the  separation  between  explosion  and  earthquake  magnitude  data.  Using  the 
regression  model,  the  Bernoulli  probability  of  an  event  being  an  explosion  is  expressed  as 
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jt(x)  = - - — —  ■  (2) 

l  +  ea+-  * 

For  n  observed  events,  earthquakes  and  explosions,  the  likelihood  function  is  defined  as  the  product  of  these 
probabilities: 


t(a,(5)  =  )*'  (1  - 

/=1 


where  y,  =  1  if  the  event  is  an  explosion,  y.  =  0  if  the  event  is  an  earthquake,  and  x,  is  the  vector  of  observed 
magnitudes  for  the  event.  Maximizing  t(cL,f5)  provides  estimates  (calibration  values)  for  a,  /3.  For  a  new  event, 

Jt(X)  is  evaluated  with  magnitudes  X  and  identification  made  with  this  value.  For  example,  if  Jt(X)  >  0.55  the 
event  is  identified  as  explosion,  if  Jt{X)  <  0.45  the  event  is  identified  as  earthquake,  and  indeterminate  otherwise. 
Figure  6a  gives  the  function  7t(X)  , 

n{EX\M ^Rayleigh),  Ms(Lovej)  =  1+e(4.09+12.14.Mj(to^,)-lz.6S.Mj(Jinytetoh),  (4) 

using  the  average  jackknife  parameter  values  (see  Figures  6b,c,d),  and  a  subset  of  jackknife  event  identifications. 
The  indeterminate  region  is  included  on  plot.  We  have  completed  a  leave-one-out  (jackknife)  cross  validation 
analysis  (Figures  6  b,c,d)  on  M/VMAXj-Rayleigh  and  Love  using  the  decision  rule  above.  The  data  included 
82  explosions  and  264  earthquakes.  For  each  jackknife  sample,  82+264=346  in  total,  calibration  values  for  Of,  /3 

were  computed  using  maximum  likelihood  estimation.  These  parameter  values  were  then  used  to  identify  the 
hold-out  event  by  evaluating  Jt(X)  and  applying  the  decision  criteria  above.  The  performance  of  the  cross 
validation  analysis  is  given  in  Table  1.  From  Figure  6e,  we  note  that  the  absolute  values  of  the  slopes  for 
earthquakes  and  explosions  are  statistically  different. 


Table  1.  Cross  Validation  Identification  Performance  with  Ms  Rayleigh  and  Love  Magnitudes. 


EX 

EQ 

I 

Total 

EX 

57 

22 

3 

82 

EQ 

11 

246 

7 

264 

The  cross-validation  analysis  of  the  proposed  MS:MS  discriminant  correctly  identifies  57  of  82  explosions  and  246  of 
264  earthquakes.  The  analysis  misidentifies  22  explosions  as  earthquakes  and  11  earthquakes  as  explosions.  These 
results  show  that  there  is  discrimination  information  in  an  M/Rayleighj.'MdLove)  discriminant.  Further  comparative 
research  is  planned  for  M/Rayleigh)  :M,(Love)  versus  Ms:mb  using  common  data.  We  fully  expect  that 
Ms(Rayleigh):Mj(Love)  will  contribute  significantly  to  multivariate  event  identification.  Results  from  this  study  do 
suggest  that  a  Ms-Love:mh  discriminant  might  be  more  robust  than  A7s-Rayleigh:m/,  due  to  the  typically  larger 
MrLove  magnitudes  for  earthquakes  and  smaller  values  for  explosions.  However,  the  smaller  M,-Love  estimates  for 
explosions,  while  great  for  discrimination,  are  costly  in  terms  of  detection.  Bonner  et  cil.  (2006)  determined  for 
MS(VMAX)-Rayleigh  to  be  measured  at  the  NNSS,  the  mb  must  be  3.6  or  greater;  thus  the  event  body-wave 
magnitude  for  M,(VMAX)-Love  application  would  increase  to  greater  than  4.0. 

Possible  Improvements  to  Love  Wave  Magnitude  Estimation 

The  objective  of  this  study  was  to  evaluate  Equation  1  as  is  for  Love  waves  from  both  earthquakes  and  explosions; 
however,  we  do  note  that  some  terms  in  Equation  1  could  change  for  Love  waves.  In  the  next  few  paragraphs,  we 
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discuss  possible  changes  for  future  application  of  this  technique,  including  the  source  excitation  and  attenuation 
corrections  and  the  need  to  incorporate  additional  periods  into  the  analysis. 


12.0  12.2  12.4  12.6  12.8  13.0  13.2  13.4 


b)  Mean=4. 08668  Variance=0.00206854 


c)  Mean=-12.646  Variance=0.00638099 


-13.4  -13.2  -13.0  -12.8  -12.6 

Mean=12.137  Variance=0.00642287 


Figure  6.  Logistic  regression  results  for  a  possible  Ms(Rayleigh):Ms(Love)  discriminant,  a)  using  the 

average  of  jackknife  parameter  values,  b)  Jackknife  O.  estimates,  c)  Jackknife  estimates  of  Rayleigh 
slopes,  d)  Jackknife  estimates  of  Love  slopes,  e)  Histograms  of  the  absolute  value  jackknife  slopes 

— .  Rayleigh  jackknife  slopes  are  gray  and  Love  jackknife  slopes  are  light  gray. 

Excitation  Correction.  The  “source  excitation”  correction— 0.66  log  f-jpj  in  Equation  1  is  slightly  misleading, 
because  the  correction  is  actually  for  the  effects  of  the  source  depth  and  structure  at  the  source  rather  than  for  the 
actual  source  spectra.  A  shallow  explosion  will  generate  large  amplitude,  short-period  (<  20  secs)  surface  waves 
(and  magnitudes)  relative  to  To=20  seconds,  where  most  historical  measurements  have  been  made,  and  thus  must  be 
reduced  in  order  to  improve  explosion  and  earthquake  discrimination  and  provide  better  agreements  with  historical 
magnitude  scales.  This  correction  accomplishes  this  need  and  is  determined  empirically  by  modeling  Rayleigh 
waves  generated  by  1  km  deep  explosions  in  a  variety  of  different  velocity  structures  (Bonner  et  ah,  2006). 

Since  Love  waves  are  not  generated  by  isotropic  explosions,  we  determined  a  corresponding  correction  for  the  Love 
waves  using  a  1  km  deep  double-couple  earthquake  for  different  velocity  structures.  A  similar  expression  for  the 

Love  wave  source  excitation  is  —0.45  log  which  is  similar  to  the  one  incorporated  into  the  Russell  (2006) 
equation.  Lor  future  examination,  a  more  rigorous  and  model-dependent  approach  to  this  correction,  such  as 
discussed  in  Stevens  and  McLaughlin  (2001)  and  Stevens  et  al.  (2007),  could  lead  to  improved  results. 
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Attenuation  Correction.  To  investigate  the  applicability  of  the  attenuation  correction  in  Equation  1  to  Love  waves, 

fT  \  1-8 

we  first  subtracted  the  correction  [0.0031  (  —  1  A]  from  our  estimated  Ms{  V VI A X )  then  computed  new  corrections 

in  the  form  [crA]  for  both  Rayleigh  and  Love  waves.  The  attenuation  coefficients  calculated  as  a  result  are  a  = 

0.0037  for  the  Rayleigh  and  a  =  0.0042  for  the  Love  waves,  compared  to  0.0031  in  the  original  formula.  Application 
of  the  new  attenuation  corrections  improves  the  residuals  for  the  events  used  in  the  inversion;  however,  it  did  not 
improve  the  RMS  residuals  for  the  entire  Middle  Eastern  data  set.  Future  application  of  this  technique  could 
possibly  incorporate  2D  or  3D  attenuation  models  for  Love  and  Rayleigh  waves  (Levshin  et  ah,  2007;  Stevens  et  ah, 
2006,  2007). 


Period  Limitations.  M/VMAX)  was  originally  designed  for  Rayleigh  waves  in  the  period  range  between  8  and  25 
sec.  However,  for  some  focal  mechanisms,  deep  events,  and  complex  paths,  the  maximum  of  the  surface  wave 
amplitudes  may  be  achieved  at  longer  periods.  Limiting  the  period  range  to  25  seconds  for  the  Middle  East  data 
(Ligure  4)  results  in  numerous  measurements  that  are  “pegged”  at  the  upper  limit  of  25  sec.  Increasing  the  upper 
limit  from  25  sec  to  40  sec  results  in  higher  magnitude  measurements  (Figure  7;  compare  to  Figure  4a),  which 
provides  a  more  reliable  estimate  of  source  size  and  slightly  lower  standard  error  for  both  Rayleigh  (e.g.,  0.22  m.u. 
to  0.21  m.u.)  and  Love  waves  (0.22  m.u.  to  0.20  m.u). 


Ms  (25)  Rayleigh 


Ms  (25)  Love 


Figure  7.  Results  of  extending  to  analysis  periods  for  MS(YMAX)  to  40  secs  for  a)  Rayleigh  and  b)  Love  waves 
in  the  Middle  East.  The  estimated  c)  Rayleigh  and  d)  Love  wave  magnitudes  are  often  increased  by 
extending  the  analysis  period  to  40  secs. 

CONCLUSIONS  AND  RECOMMENDATIONS 


We  conclude  that  estimating  a  Love  wave  magnitude,  using  the  same  formula  and  methods  employed  for  Rayleigh 
waves,  can  lead  to  improved  earthquake  and  explosion  discrimination  due  to  the  fact  that  the  earthquakes  typically 
have  a  larger  M,-Love,  while  explosions  normally  exhibit  a  smaller  MrLove  when  compared  to  the  Afs- Rayleigh. 

We  conclude  that  an  MS:MS  discriminant  is  possible;  however,  it  does  not  have  the  same  population  separation  that 
has  been  historically  observed  for  the  Ms:mb  discriminant.  Results  also  suggest  that  incorporation  of  Love  waves  into 
the  analysis  requires  a  re-examination  of  the  period  limits  currently  used  for  the  M/VMAX)  technique.  While  the 
M s(  VM  AX)  -  Ray  leigh  method  is  currently  operational  at  different  data  centers  using  periods  between  8  and  25 
seconds,  we  do  believe  that  future  processing  should  be  extended  to  40  seconds,  especially  in  regions  with  deep 
earthquakes  and  complex  paths. 
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